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1 Introduction 

Despite the tremendous success of quantum chromodynamics (QCD) in describing the phenomenol- 
ogy of strong interactions of elementary particles at high energies, a deep understanding of the physics 
at low energies is still lacking. At high energies, non-Abelian gauge theory is asymptotically free, and 
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as a result, perturbation theory is an adequate tool. However, at low energies, the interaction be- 
comes strong, and perturbation theory fails. It has been suggested that this failure is directly linked 
to the phenomenon of confinement in QCD, which is a well-known experimental fact. However, the 
precise nature of a non-pertubative mechanism ensuring confinement is still not well understood. In 
field-theoretic terms, this means that the vacuum of QCD at low energies has a far more complicated 
structure than the trivial perturbative one. 

A model of a non-perturbative vacuum for an SU{2) gauge theory was put forward in 1977 by 
Savvidy [16]. He proposed an explicit ansatz for the vacuum gauge fields in form of a constant 
chromomagnetic field, or more precisely, a gauge field with covariantly constant field strength in 
flat four-dimensional Minkowski space-time with only one nonvanishing color component. Savvidy 
showed that, due to quantum fluctuations of the gauge fields, the energy of such a field configuration 
is below the perturbative vacuum level, which leads to infrared instability of the perturbative vacuum 
under creation of a constant chromomagnetic field. Further investigations [13, 14] showed that the 
Savvidy vacuum itself is unstable too, meaning that the physical nonperturbative vacuum has an even 
more complicated structure. It has been suggested that the real vacuum is likely to have a small 
domain structure with random constant chromomagnetic fields (spaghetti vacuum). 

In our papers [3, 4] we extended Savvidy's investigation by considering more complicated gauge 
groups and flat spacetimes of dimension higher than four. We showed that for an arbitrary compact 
simple gauge group in dimensions higher than four there exist more general nontrivial field configu- 
rations with several color and space-time components that turn out to be stable. In [4] we proposed 
an explicit example of such background field configurations. 

In the present paper we propose a new mechanism to stabilize the Savvidy vacuum in four di- 
mensions. The main idea of this approach is that a positive space curvature could provide an effective 
mass term for the gauge fields on the chromomagnetic vacuum, thus, making the vacuum stable. To 
simplify the calculations, we consider the space-times with compact space slices with the product 
structure x S^. To also study the finite-temperature effects we consider Euchdean spacetimes of 
the formal xS^xS^. 

However, as we show below, topological considerations on the sphere constrain the magnetic 
field to be of the same order of magnitude as the space curvature, thus negating the stabilization effect 
of the curvature term. Moreover, even under these constraints, an interesting second-order phase 
transition occurs at a critical temperature near the inverse radius of the sphere. 

2 Yang-IVIills Theory 

Let (M, g) be an ^-dimensional pseudo-Riemannian orientable spin manifold without boundary with a 
globally hyperbolic metric g. We denote the local coordinates on M by ;c^, with Greek indices running 
over 0, 1, 1. We denote the frame indices by the low case Latin indices from the beginning 
of the alphabet, which also run over 0, 1, . . . , « - 1. The frame indices should not be confused with 
the group indices introduced below that are enclosed in parenthesis. We use Einstein summation 
convention and sum over repeated indices. The coordinate indices are raised and lowered by the 
metric tensor gjjy and the frame indices are raised and lowered by the Minkowski metric, rjat- 
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We choose a local Lorentz frame on the tangent bundle TM, Ca - ^a^d^, and the dual frame on 
the cotangent bundle T*M, cr" = cr"^dx'^. We denote the corresponding spin connection 1-form by 
aj"b - oj"bfidx'^ and its curvature 2-form by 0";, - ^©"h^ydx^^ A dx^ so that Rev - ea'^e'^a®"bnv is 
the Ricci tensor, and R = g'^'^Rfjv is the scalar curvature. Let 7~ be a spin-tensor bundle realizing a 
representation T of the spin group, Spin(l,« - 1), with generators Xah- The spin connection induces 
a connection on the bundle T with the curvature "Rfjy = ^&"''fjvT(Lab) ■ Recall that the generators 
of the spin group in the vector representation and the spinor representation are 

TlC^abfd - '26'^[arib]d , T'spinCSflfo) - - Jab, (2-1) 

where t/m is the Minkowski metric and yah - yiajhy Here and everywhere below the square brackets 
denote the antisymmetrization over indices included. 

Let G be a m-dimensional compact simple Lie group and g be its Lie algebra. We use lower case 
gothic letters to denote Lie algebras, for example, spin(l,« - 1). We denote the group indices, which 
run over \,2, . . . ,m„ by the low case Latin letters from the middle of the alphabet. Let C ju be the 
structure constants of G in a given basis. The real m x m matrices C, defined by (C,)\. = C^tk form 
a basis in the Lie algebra g and define the adjoint representation ad : g ^ End(R'") of the algebra 
g by endomorphisms of the vector space R'". Of course, this also defines the adjoint representation 
Ad : G — > Aut(]R'") of the group G into the automorphism group of the vector space R"'. In the 
following we will identify the algebra g with its adjoint representation. The Cartan-Killing metric jij 
on the Lie algebra g is defined by 

irCiCj = C\,djk = -2yij, (2.2) 

the precise form depends, of course, on the choice of the basis. We will determine it later. For compact 
semi-simple groups it is positive definite. We will use it to raise and lower the group indices. 

We consider the principal fiber bundle over the manifold M with the structure group G and the 
typical fiber G. Let p-w : G — > Aut(VK) be an irreducible representation of the group G into the 
automorphism group of an A'^-dimensional (real or complex) vector space W. Sometimes, we will 
denote the representation p-w by the vector space W; this should not cause any confusion. For example, 
we will denote the generators of this representation by W{Ci). Let be the associated vector bundle 
with the structure group G and the typical fiber W. Then for any spin-tensor bundle T (realizing a 
representation T of the spin group, Spin(l, n - 1)) the vector bundle "V = "W (g) 7" is a twisted spin- 
tensor bundle realizing the representation V = W ®T . The sections of the bundle ai^e represented 
locally by (real or complex) A'^-tuples of spin-tensors. 

Let = A'^Cj and ^ = ^^dxf^ be the Yang-Mills connection 1-form taking values in the Lie 
algebra g and 7^,y = T'^yCi and T - ^ffiv dx^ A dx^ - dJ{ -i- A ^ be its curvature 2-form. It is 
worth remembering that in a non-trivial bundle, there ai^e several overlapping coordinate patches; the 
connection one-form is not globally well-defined, in general, but is rather described by a collection of 
its representations in each patch, which are related in overlapping patches by gauge transformations. 

Let be the total connection on the twisted spin-tensor bundle "V. The curvature of the total 
connection on the twisted spin-tensor bundle is equal to ly/®T{'Rf^y)+W{'F'ny)®lT . We will usually omit 
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the identity matrices for the sake of simphcity of notation. The covariant Laplacian t^r^w - g^^^^^X 
acting on sections of the twisted spin-tensor bundle has the form 



AtkiW - g 



dy + -co' 



cd 



Tilled) + 



(2.3) 



We will use twisted Lie derivatives defined as follows. Suppose that there is a faithful represen- 
tation px ■ Spin(l,?i - 1) ^ Aut(VK) of the spin group in the same vector space W with generators 
X(Lab)- Then the matrices 

Gab = V ® TiZah) - XiZah) ® Ir (2.4) 

are the generators of the twisted representation pw»T '■ Spin(l,?i - 1) ^ Aut(V) of the spin group. 
Let ^ be a Killing vector field. The twisted Lie derivative of sections of the vector bundle 'V along ^ 
is defined by 



f, 2^[a;b]G' 



ab 



= + \ [eojah^ - ^[a-M] n^"') + eA'^^wid) + y^a^xii:"'') . (2.5) 

The action of the Yang-Mills theory in curved spacetime is constructed as follows. We consider 
two associated vector bundles "IVspin and "Wq of dimensions A^spin and A^o respectively realizing some 
irreducible representations, Wspm and Wq, of the gauge group. Usually, the representation Wspin re- 
alized by the spinor fields is taken to be the fundamental (or defining) representation of the gauge 
group. The scalar fields are just sections of the bundle "Wq, whereas the spinor fields are sections of 
the spinor bundle twisted by "IVspin- Then the classical action of the model is the functional 



M 



(2.6) 



where g - det^^y, \T\^ - AdT^vF'^^ , e is the Yang-Mills coupling constant, Vitp) is a potential 
for scalar- fields (such that V(0) = V(0) = 0) and M is the spinor mass matrix. 



3 Effective Action 



In this paper we will not quantize gravity, assuming the gravitational field to be classical and simply 
ignoring all quantum-gravitational effects. The energy scale of our primary interest will be well below 
the Planckian scale, so that this assumption is reasonable in any theory of quantum gravity. 

The effective action is expressed in terms of functional determinants of differential operators 
(see, for example, [9, 10]). Such determinants can be defined, strictly speaking, only for elliptic 
operators on compact manifolds by making use of a regularization procedure, for example, zeta- 
function regularization. For hyperbolic operators and for non-compact manifolds such determinants 
do not have a direct mathematical meaning. 
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That is why we will make the further assumption that the background is static, that is, there is a 
global time-like Killing vector field df Moreover, we will assume that the spacetime has the simple 
structure of a product manifold M -M.x'L and that all background fields are static and do not have 
time-like components. 

In the case of static background one can make an analytic continuation to a purely imaginary 
time t h , with a positive-definite Riemannian metric. Moreover, we can go even further and 
compactify the Euclidean time by replacing R by a circle of radius a\, that is, by restricting the 
range, < t < yS, where p = Inai is the circumference of the circle S ' and requiring all fields to be 
periodic in the Euclidean time r with period /?. The "Euclidean" space-time M = 5 ' x Z is then a 
compact manifold and Lorentzian spin group Spin(l, «- 1) becomes the Euclidean spin group Spin(?i), 
which is compact. This corresponds to a statistical ensemble at a finite temperature T = 1 //?. In the 
limit of infinite radius jS ^ oo we recover the zero-temperature theory. 

We will consider a background in which there are no matter fields and Yang-Mills fields and 
gravitational field are covariantly constant (parallel), that is, 

V^Vt»/3 = 0' V^r„/3-0. (3.1) 

More precisely, we will study the case when E = 5 ' x S^, that is, M = x x S^. 

There exists a minimal gauge such that all differential operators involved are second-order differ- 
ential operators of Laplace type 

L = -A + Q (3.2) 

with some endomorphism Q. The precise nature of the Laplacian A = g'^^V^Vy depends, of course, 
on the fields it is acting upon. 

Let L^^^^ be the Laplacian acting on scalar fields in adjoint representation of the gauge group 
defined by 

-^gtost ^ -^To^Ad ■ (3.3) 

Let Qvect be the endomorphism acting on vectors in adjoint representation of the gauge group defined 
by 

(Gvect^)" - (RMAd - 2r\) if'' . (3.4) 

Let Lvect be an operator acting on vector fields in adjoint representation of the gauge group defined by 

■^-vect - -^Ti^Ad + Gvect ■ (3.5) 

We will suppose that the spinor mass matrix M commutes with the Dirac operator 7^ V^. Let 2spin 
be the endomorphism acting on spinor fields in the representation Wspin of the gauge group defined by 

Gspin - \rIt,,,„ ® H.„ - Ir^W^p^niTat) + . (3.6) 

Let Lspin be the differential operator acting on spinor fields in the representation Wspin of the gauge 
group defined by 

^spin = -^r^pin^Wspin + Gspin • (3.7) 
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Let Lo be a differential operator acting on scalar fields in some representation Wq of the gauge 
group defined by 

Lo - -/^To»Wo + Qo ' (3-8) 



where Qo is a matrix defined by 



„ - {h, Qoh)wo ■ (3-9) 

e=0 



The most important observation that should be made at this point is that the positive curvature 
acts as a mass (or positive potential) term in both the Yang-Mills operator and the spinor operator. 
While the magnetic field reduces the eigenvalues of the Yang-Mills operator the positive Ricci tensor 
increases them. Roughly speaking, it is the balance of these two terms that determines whether or not 
the Yang-Mills operator is positive (so that the vacuum is stable). 

In the minimal gauge the one-loop effective action is given by [3, 4] 

Y = S+h (r(i)yM + r(i)n,at) + 0{h^) , (3.10) 

where 

T(\)YM = 2 log Det Lvect - log Det ^ , (3.11) 
and ^ ^ 

r(l)mat = 2 log ^° ~ 2 ^'P'" ' 

axe the contributions of the Yang-Mills field, and matter fields respectively, and Det is the functional 
determinant. The Planck constant is introduced here just for illustrative purposes. Henceforth, we set 
h=l. 

The operators introduced above are second-order elliptic partial differential operators on compact 
manifold. The spectrum of these operators depends, of course, on the background fields. Elliptic 
operators on compact manifolds can only have a finite number of negative eigenvalues. The negative 
eigenvalues indicate instability of the vacuum at low energies. 

That is why, to study the infrared behavior of the system, one has to introduce an infrared reg- 
ularization and take it off at the very end as in [4]. For example, one could introduce a sufficiently 
large mass parameter z so that all operators are positive, which is equivalent to replacing the operators 
Lhy L + z, and study the dependence of the infrared regularized effective action on z- If there are no 
infrared divergences then there is a well defined limit z ^ 0. In the case of non-trivial low-energy 
behavior there appears an imaginary part of the effective action or some infrared logarithmic singu- 
larities. One should stress that although the ultraviolet regularization is a rather formal method, the 
infrared regularization parameter can take, in principle, a direct physical meaning, something like 

We will assume that this has been done, that is, we will add a mass parameter to the Yang- 
Mills and the ghost operators so that all operators are positive. The determinants of positive elliptic 
operators can be regularized by the zeta-function regularization method which can be summarized by 
[1,5,9] 

logDet(L + z)--r(0), (3.13) 
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where ^(0) - IMs) 



and 

.v=0 _ 



where = AkiL) are the eigenvalues of the operator L and dt are their multiplicities. The analytic 
continuation of the zeta function gives a meromorphic function of s, which is analytic at s = 0; 
therefore, the determinant (3.13) is a well defined invariant. 

We can also express the zeta-function in terms of the heat trace of the operator L defined by 

oo _ 

Tr exp(-fL) = V dkc"''^' - I dx g^'hr uf^^{t) , (3.15) 

where Uiit', x, x') = exp{-tL)6{x, x') is the heat kernel of the operator L and U'j^'^^if, x) - Ui{t; x, x) 
is the heat kernel diagonal. Then the zeta-function is related to the heat trace by the Mellin transform 



^{s) Jdt f'^e~''Tr exp(-?L) , (3.16) 





where /lis a. renormalization parameter introduced to preserve dimensions. 
We introduce a useful function 0z,(O as follows 

0l(O - {4nt)"'^Tr exp(-fL) . (3. 17) 

It is well-known that as ? ^ 

oo 

0z.(O~2]s/, (3.18) 

k=0 

where = Bk{L) are some spectral invariants of the operator L. The first three coefficients have the 
form [1,5, 10] 

Bq= \ dxg^'hr:!, (3.19) 
Jm 

Bi = dxg^'hr: {^Rl - , (3.20) 

+ ^VRab+Tah^[^"'' + T"'''^. (3.21) 

It will be convenient to represent the heat trace as ? ^ as follows 

oo 



*=0 
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where i is a new arbitrary parameter (that should not be confused with z) and 

k 



7=0 



(3.23) 



The analytic continuation of the zeta function to 5 = can be obtained by integration by parts. 
Then the zeta-regularized determinant can be expressed directly in terms of an integral of the heat 
trace. For even n we obtain 



log Det^(L + z) - [logOu'O + ^ (1 + ^)| (^) l^~"®Lit)] , (3.24) 



where ^(s) = T'{s)/T{s) is the logarithmic derivative of the Gamma function. 
Let us consider the case of four dimensions, « = 4, in more detail. Then 



logDet^(L + z) 



dt 



1 3 
logOu^O + - - I 



1 1 [e-"&L{t)\ , (3.25) 



dt 



where C ^ 0.58... is the Euler constant. It is not difficult to find the dependence of the determinant on 
the renormalization parameter ju 



logDet/,(L + z) - -(47r)-2log^ A2(-z) + logDet ^{L + z) . 

Now, let 

OrCO = - 1^0 + Ay{A)t + A2{A)t^] , 
The renormalized heat trace has the following asymptotics: as f — > 

&l\t) = A3{A)t' + 0{t^) , 

and as ? ^ 00 

&l^(t) = -e"'^ [A2{A)f- + Ai{A)t + Ao(i)] + {Anf dit^ e''^' + 0{e''^^) . 
Thus, we can consider the integral 

Jr*CO 1. 
-e~''&l\t), 
f 



(3.26) 



(3.27) 



(3.28) 



(3.29) 



(3.30) 



which is well defined since it converges both at and 00. 

One can compute the dependence of the renormalized determinant on A exactly. First, we show 

that 

log Det,en(Z. + z) = -{AnT^ Uii-z) + AA,{-z) + . (3.31) 
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Then, by integrating this equation we get 

^ . . , , 1 



logDetren(Z. + Z) - -{471)' 



Aq 4 



+ const , 

(3.32) 



where Aq is some constant. Notice that in the hmit when A ^ there is an infrared divergence 

logDet,en(^ + z) - -(47r)-2A2(-z)log ^ + 0(1) . (3.33) 
Now, by integrating by parts one can show that 

logDet^(L + z) = logDetren(Z. + z)-(47r)-2log^ A2(-z) + coi2Ao + ci^i(-z) + C2A2(-z). (3.34) 

Here co,c\ and C2 ai^e some numerical constants dependent on the regularization scheme, in particular, 
they can be set to zero without loss of generality. 

By using this regularization of functional determinants we obtain the effective action in the form 



^WYM = -2(4^)"' {/^YM log J + l J e~'Wy'Ut) \ , (3.35) 



where 



^(D^at = -^(4^)"' {>Smat log ^ + J^^ f ©matCoj , (3.36) 

0™(O = 0ri(O - (0 , (3.37) 

ghost 

0Z(O = 0^"(O-0r„(O, (3.38) 



and 



Z2 

J3yM - 52(Lvect) - zBi{Lyect) + -^BoiL^sct) 

-2B2{L™^;) + 2zBi(Lghost) - z'Bo(ighost) , (3.39) 

)8mat - 52(Lo) - B2(Z.spin) - (3.40) 

The main idea of the renormalization group is based on the realization that the total effective 
action, T = S + TiT^i) + • • ■ , should not depend on the arbitrary renormalization parameter /u. This 
means that in renormalizable field theories the coupling constants in the classical action should depend 
on yu in such a way to exactly compensate the dependence of the one-loop effective action on /i, that 
is, 

li^S = Fd) = (4;r)-2 (J3ym + Anat) ■ (3.41) 
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This means that the classical action should have terms of the same type as those in the coefficients 
Pym and/3niat- 

As we discussed above, all relevant operators are of Laplace type L = -A + Q. For a covariantly 
constant background the endomorphism Q is covariantly constant, and, therefore, commutes with the 
Laplacian. Therefore, the heat semigroup of the operator L is determined by the heat semigroup of 
the Laplacian 

exp(-?L) = exp(-?2) exp(?A) , (3.42) 
so, the heat kernel diagonal of the operator L has the form 

uf'Ht) = t^vi-tQW'^^Kt) , (3.43) 

where ?7'^'^§(?) denotes the heat kernel diagonal of the pure Laplacian acting on a vector bundle "V. 

By using this property we can express the heat kernel diagonals of the operators introduced above 
in terms of the heat kernel of the corresponding Laplacians, 

^ty^ = ^'^P (-^2vect) , (3.44) 

t^?M(0 = C(0, (3.45) 

ghost 

Ul'\t) = exp i-tQo) t/?i;V„(0 , (3.46) 
<!« = -P(-^M<U,„«- (3.47) 

This means that 

@YMit) = vol (M)(4;rO'tr^4tr r, exp(-tQ,,,OUf^^^/t) - 2[/J5^/o] , (3.48) 
0mat(O - vol(M)(47rO'|trwoexp(-feo)t/JX(0 

^diag 

spin'^^'spin 



-tr w.„„tr c^vi-tQ,,,,)U-^2^ St) \ . (3.49) 



The renormalized functions &y^(t) and ©]^"t(^) obtained from this by finding the smallest eigen- 
value and then subtracting some terms according to the prescription (3.27). Thus, we need to compute 
the heat kernel diagonals for the Laplacians only. 



4 Geometry of the Sphere S ^ 

In this section we follow mainly our paper [7]. 
4.1 Metric 

We cover the sphere 5^ (of radius a) by two coordinate patches: one patch covering the South pole 
and another patch covering the North pole. We will use the spherical coordinates (r, (p), which range 
over < r < an and < ^ < In. The South coordinate patch is the neighborhood of the South pole 
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r = 0, whereas the North coordinate patch is the neighborhood of the North pole r = an. The volume 
of 5^ is, of course, vol (5^) - 4na^ . 

The metric in spherical coordinates is 

ds^ - dr^ + (? ^\x?{rjd)dip^ . (4.1) 

We choose an orthonormal basis of 1 -forms 

cr' = cos ^(ir - a sin(r/fl') sin ^(i^ , (4.2) 
cr^ - sirup dr + a sm{r/ a) cos (fi dip . (4.3) 

Then the spin connection one-form is 

oJab = SabU - cos(r/a)] d(p, (4.4) 

the Riemann curvature is ^ ^ 

Rabcd - -^£ab£cd - -^{^ac^bd " ^ad^bc) , (4.5) 

and the Ricci tensor and the scalar curvature are 

Rab - -^Sab , R=\- (4-6) 

4.2 Connection 

In two dimensions the Yang-Mills connection with a covariantly constant curvature is necessarily 
Abelian and can always be chosen to be proportional to the spin connection, that is, 

J{ = -X{\- cos(r/a)] dip , (4.7) 

where X = ad(E) is an m x m real anti-symmetric matrix. The Yang-Mills curvature 2-form is 

T ^-X- sm{r/a) dr hdip, (4.8) 
a 

in components, T^d = ^cdH , where H = -X/a^ . The connection has to be redefined (via a gauge 
transformation) in the North coordinate patch. By defining 

=X[l+cos{r/a)]dip, (4.9) 

we obtain the connection that is well defined globally. Now, we have that 

^' -Jl^dUU-\ (4.10) 

where U = exp (IXip) is an element of the group G. 

Since U should be periodic in ip, the matrix X should satisfy the condition 

exp(47rX) = I, (4.11) 
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which means that its eigenvalues must be either zero or imaginary half-integers. Only in this case the 
vector bundle is globally defined. 

Since X is anti-symmetric, its non-zero eigenvalues must appear in pairs. That is, the spectrum 
of the matrix X must be 

0^^, /y, -/y , . . . , i^, -i^ , (4.12) 

where r = m - 2p, and nj, j = I, . . . , p, are some non-zero integers. 

Such construction for each n is nothing but the Hopf complex line bundle 'Hn over S^, which is 
equal to a tensor product of n Hopf bundles 'Hi (or dual Hopf bundles for n < 0). This is sometimes 
called the topological quantization condition. The deep fundamental reason for this condition is the 
Chern theorem [11], which in particular says that the Chern form 



271 Is 



T = -2iX (4.13) 

has integer eigenvalues. Note that the matrix X is the generator of a representation px : Spin(2) — > 
Aut(R"') of the double cover of the group S 0(2) (which we, by definition, still call the spin group 
Spin(2) in two dimensions). 

It is important to understand that the numbers nj are not independent. Because X lies in the 
adjoint representation of the compact semi-simple Lie algebra g, these eigenvalues are determined by 
the roots of the algebra g. We will discuss this in the next section. 

4.3 Roots 

In this subsection we follow [12]. Since X is a real antisymmetric matrix it can be diagonalized, 
and thus lies in the Cartan subalgebra. Let r be the rank of the group G, which is equal to the 
dimension of the Cartan subalgebra. The generators Cj, j = 1, . . . , r, of the Cartan subalgebra in 
adjoint representation are m x m diagonal matrices of the form 



Cj = diag 



0, . . . , 0, iaij, -iaij, . . . , icupj, -io:pj, 



(4.14) 



where at = {auj), k - 1, . . . , p, are some covectors in W called the positive roots of the Lie algebra g, 
p is the number of positive roots related to the rank by r = m -2p. Then the matrix X has the form 



X 



^ ^ x^Cj = diag 



i=i 



0, . . . , 0, iai{x), -iai{x), . . . , iap{x), -iap(x), 



(4.15) 



where x = (x^), j = 1, . . . , r, is a vector in W, and ak{x) = ^[^j a^jX^ is the canonical value of the 
covector ak on the vector x in W. The vector x must be such that the value of each root on it is a 
half-integer, that is, 

ak{x) = —, nteZ, k = \,...,p. (4.16) 
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Thus, the number of (possibly) non-zero eigenvalues is equal to the number p of positive roots of the 
algebra, and the eigenvalues themselves are equal to the values of the roots on the vector x. 

The roots aj^ and the vector x are vectors in W. However, it is convenient to consider this space 
W as the hyperplane in R'"^^ orthogonal to the vector b - ? where Cj, i - I, . . . ,r + 1, is 

the canonical orthonormal basis in W^^^. Then the roots ak and the vector x can be represented in 
terms of the basis e, of W^^ as ;c = Y.'Hi ^'^i ' ^^id = T^'i^l Q'fa-^/ , k = I,. . . ,p . Notice that not 
all of the coordinates x' are independent. Since these vectors he in the hyperplane orthogonal to the 
vector b, the sum of all coordinates of these vectors should be equal to zero, that is, f -0,and 
T^'i^l a^ki = . Then the values of the roots on the vector x are 

r+l r 

akix) ^ ^ akix' = ^ i&ki - ak,r+i) x' . (4.17) 

For all classical Lie algebras A„, S„, C,, and D„ the coordinates of the roots are integers (up to a 
uniform normalization factor), that is, = A.^ki , k = I, . . ., p; i - 1, . . . ,r + 1 , where Pki are 
integers and /I is a normalization factor. The normalization constant can be determined from the 
chosen Cartan-Killing metric (2.2) by requiring 

p 

A^Y,l^kiPkj = yij. (4.18) 
k=\ 

This can also be viewed as the definition of the metric jij. The constant X is then just a uniform factor 
that can be set to /I = 1 . 

Therefore, ak{x) will be half-integer if x' - jki , where ki,i - 1, . . . , r, are arbitrary integers, that 

is, 

ak{x) = - 2 {J3ki - A.r+i) ki . (4. 19) 

^ ,=1 

To be specific let us consider the group G - SU{N). The algebra su(A'^) is isomorphic to the 
classical algebra Aat-i. The dimension and the rank of SU{N) are 

m = dimSU{N) ^ - I, r - rankS UiN) = N-l. (4.20) 

The positive roots are labeled by two integers I < i < j < N, and have the form ajj = et - ej . The 
number of positive roots is 

m-r NjN - 1) 
P- — = ^^- (4-21) 

The Cartan-Killing metric is yij = 25 ij - - Sij+i . Then 



aijix) = x' - xJ = ^{ki - kj) . (4.22) 



Therefore, in this case the integers determining the eigenvalues of the matrix X are also labeled 
by two indices 

riij = ki-kj, I <i < j <N. (4.23) 
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We will call these integers monopole numbers for njj > (or antimonopole numbers for njj < 0). Here 
ki,i = I,. . . ,N, are ai^biti'ary integers whose sum is equal to zero. In particular-, some or all of them can 
be equal to zero. Note, however, that it is impossible to have only one non-zero number kj. Therefore, 
either they are all equal to zero, or there are at least two non-zero integers Another important 
observation is that for any choice of non-zero integers ki some of the integers «,y will have absolute 
value greater or equal to 2. In other words, it is impossible to have = 0, ±1 for all /, j (except, of 
course, the trivial case when all integers ki = 0). This observation has profound implications for the 
stability of the chromomagnetic vacuum studied in this paper. 

For the classical groups the situation is similar. Let, as above, kj, i - 1, ... ,A'^, be an arbitrary 
collection of A'^ integers whose sum is equal to zero. Then for the group D^^y the possible monopole 
numbers are 

riij - ±ki ± kj . (4.24) 
For the group B^^i there are two possible combinations 

Uij = ±ki ± kj and - ±ki , (4.25) 

and for the group Cjv-i the possible combinations are 

Hij - ±ki ± kj and = ±2ki . (4.26) 

It is not difficult to see that in all these cases there is no choice of non-zero integers kt such that the 
only monopole numbers are 0, ±1. There will be necessarily monopole numbers with absolute value 
greater or equal to 2. 

Thus, for any of the compact simple classical groups if the matrix X has at least one non-zero 
eigenvalue, then it will have at least one eigenvalue with absolute value greater or equal to 2. 

4.4 Weights 

Now let us consider an irreducible representation g End( W) of the Lie algebra g in a A'^-dimensional 
complex vector space W. The generators of the Cartan subalgebra in this representation, W(C,), 
/ = 1, . . . , r, are A'^ X A'^ complex diagonal matrices of the form 

W{Cj) = diag{ivij,...,ivNj) , (4.27) 

where = (vkj), k = I, ... ,N, are some covectors in W called the weights of the representation W. 
Contrary to roots, the weights can be degenerate, that is, have multiplicity greater than 1, and be equal 
to zero with some multiplicity too. Then the matrix W{X) has the form 

r 

WiX) = 2 x^WiCj) = diag (/vi (x),..., ivNix)) , (4.28) 

!=1 

where x = (x^), j = 1, . . . , r, is a vector in R'', and Vk(x) = Y,'i=i ^kjX^ is the canonical value of the 
CO vector ak on the vector x in R''. The vector x must be such that the value of each weight on it is a 
half-integer, that is, 

Vk{x) = Y^ mt^Z, k=l,...,N. (4.29) 
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The weights Ue in the same space as the roots. So, we can represent them by - T,'i=i ^ki^i , 
k - I, . . . ,N , where T^'ill ^ki - ■ Then the values of the weights on the vector x are 

r+l r r 

Vkix) ^ ^ Vkix' ^ ^ (vki - V/t,r+i) ^' ^ 2^ ~ ^'^•'■+1) ^' ' (430) 

i=l i=l (=1 

where ki, i = 1, . . . , r, are arbitrary integers. 

To be specific let us consider the fundamental (defining) representation of the algebra g = su{N) 
by N xN complex traceless anti-Hermitian matrices. The Cartan subalgebra is generated by diagonal 
matrices. Therefore, the generators W{Ci), i = I,. . .,N - 1, of the Cartan subalgebra must have N 
imaginary eigenvalues whose sum is equal to zero. Then the matrix W{X) must have the form 

W{X) = g xJWiCj) = diag • ■ • , /y) , (4.31) 

where kj, i = I, . . . ,N, me N integers whose sum is equal to zero. These are exactly the integers that 
define the values of the roots aij{x) = (ki - kj)/2 for the adjoint representation of the group SU(N). 

4.5 Isometrics 

It is well-known that = S 0(3)/ S 0(2), so that S0{3) is the isometry group and S0{2) is the 
isotropy (or holonomy) group of S^. The Killing vectors, ^a, of have the form [7] 

^1 - cos (fdr- - cot(r/a) sin (f , (4.32) 
a 

^2 - sin (fidr + - cot(r/a) cos (p 5„ , (4.33) 
a 

(4.34) 

One can check that the Killing vector fields form a representation of the isometry algebra, S 0(3), [7] 

[^1,^2] = -4^3, [ft,^i] = -6 [6,6] = -^1- (4.35) 

a^ 

The Cartan metric of the group S 0(3) has the form 

(tab)- diag (l,l,fl2), (/«)-diag|l,l,^J . (4.36) 

Therefore, the Casimir operator, A = j^^^a^b, of the Lie algebra of the group S 0(3) is nothing but 
the scalar Laplacian 

A = dj + - cot{r/a)dr + ^ \ ^ dl . (4.37) 
a a^sm'^ir/a) 

Now, suppose the group G has the spin group Spin(2) as a subgroup and let a : Spin(2) — > G be 
the con^esponding embedding. Since we also have a representation pw '■ G Aut(VK) of the gauge 
group G in the vector space W, this defines a new representation of the spin group px = Pw ° or : 
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Spin(2) — > Aut(lV). Let S be the generator of the spin group Spin(2). Then is the generator of 
the spin group Spin(2) in the representation px and 

G = lw® T{Y) - ® It (4.38) 

is the generator of the twisted representation X®T of the spin group Spin(2). This generator should 
not be confused with the gauge group denoted by the same symbol. 

The twisted Lie derivatives Xa = -C^a along Killing vectors of sections of the vector bundle 
ai-e [7] 

1 1 - cos(r/a) 

Xi = cos (po,- sin w cot(r/a) do, + sin w G , (4.39) 

a a sm(r/a) 

1 1 - cos(r/a) 

X2 - sirup Or -\ — cos ^cot(r/fl') Ou, - cos ^ ■ G, (4.40) 

a asm(r/fl) 

£3^d^ + G. (4.41) 
One can show that these operators form a representation of the isometry algebra so(3) [7] 

Ul,£2] = -\£3, [-C3,-Cl] = -£2 [X2,X3]--Xi. (4.42) 

The generalized Laplacian is expressed in terms of the Casimir operators of the isometry group 5 0(3) 
and the holonomy group S 0(2) 

A = j^^'ZaZb - \g^ (4.43) 

and is equal to [7] 

A = dl + - cot{r/a)dr + \ (d^ + [1 - cos(r/fl)]G)^ . (4.44) 

a sin^(r/a) ^ ^ 

We will be interested in the limit as a ^ 00 (when the curvature of the sphere vanishes, so, 
formally, M?) and X — > 00, so that the magnetic field H - -X{l,)/a^ remains constant. In 

this limit the generator G becomes G — > a^H (g) I7-. It is worth noting that this limit is only defined 
locally in the South coordinate patch. Therefore, the global bundle structure is destroyed in this limit. 
Thus, in the following we will be working in the South coordinate patch (for small r) in S^, which 
approaches M? as a ^ 00. 

Notice that as a ^ 00 the Killing vectors of the sphere 5^ become the Killing vectors of the 
Euclidean space and the Laplacian becomes exactly the Laplacian with a constant magnetic field 
H in the Euclidean plane M?. 



5 Heat Traces 

We will employ the algebraic methods for calculation of the heat kernel developed in [6-8]. To 
compute the heat trace directly we need to compute the spectrum of the Laplacian, its eigenvalues 
and their multiplicities. However, we can also compute the heat trace as an integral of the heat 
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kernel diagonal. Since on the sphere the heat kernel diagonal is constant the heat kernel is simply 
proportional to the fiber trace of the heat kernel diagonal, 

Tr exp(-fL) - vol (M)tr exp(-?2)[/''^^§(0 . (5. 1) 

That is why, we need to compute the heat kernel diagonal of the Laplacian. 

Next, we note that on the product manifold 5 ^ x 5 ' x the Laplacian splits naturally 

^Six5'xS2 = ^S'x5i + > (5-2) 

and, therefore, the heat kernel factorizes 

^:^s.x.^«-^''^..«^?«- (5.3) 

Therefore, 

Tr exp(-fL) - vol(M)tr exp{-tQ)uf^^^,{t)uf^\t) . (5.4) 
That is why, we need to compute the heat kernel diagonals of the Laplacian on S ' and S ^ only. 

5.1 Heat Kernel on B? and xS^ without Magnetic Field 

The Laplacian on 1. is A = 5^ . The heat kernel diagonal of such an operator is easily computed by 
Fourier transform 

.diag, - 1 



For R we obviously have the product 



U-\t) = —= . (5.5) 



= ^. (5.6) 



On the circle S ^ of radius a the heat kernel can be computed by Fourier expansion. The spectrum 
of the operator (-A) is 

^/(A) - -J , (5.7) 

where / = 0, 1,2, ... , with multiplicities do = I and di = 2 for / = 1,2, ... . Then the heat kernel 
diagonal on S Ms 

t/,r(0 = ^n(4), (5.8) 



where 



and 6T,{v,q) is the third Jacobi theta function. Therefore, for the product 5 ' x S ' with radii a\, a2 the 
heat kernel diagonal is 



^six5'('^ 47rf 



1 







' t ' 




a 















(5.10) 
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An important property of Q{t) is the Poisson duality formula, which gives a nontrivial relation 

(5.11) 



By differentiating the equation (5.1 1) we get 

2;r2 



(5.12) 



which immediately shows that Q is an increasing function. 

We will need the asymptotics of this function. It is not difficult to see that as f ^ 



0(0 = 1 +2exp|-yJ + 0{e-^'''/') , 
a'Xt) = 2nh-^ (n^ - 2t) exp (-7) + ^ (e-^""'") , 



and as f ^ 00, we have 






0(0 = 


J- 


■^1/2 


O'(0 = 


1 

-\/n 


1 

r 


O"(0 = 


1 

■y/n 


1 

-4' 



L~'l'-^(t-'l^-2t'l^)e-' + 0(e-'') 



-3/2 



+ \2t'^^-2t 



'1/2 



(5.13) 
(5.14) 
(5.15) 

(5.16) 
(5.17) 
(5.18) 



Notice that although, in general, as f ^ 00 the derivatives of the function 0(0 decrease only as 
powers of t, a particular combination of the derivatives, which we will need later, is exponentially 
decreasing as f ^ 00, 



O'(0 + 2fO"(0 = —t^'^ {2t - 3)e'' + 0(e"^') . 



(5.19) 

This will have a major impact on the calculation of the heat capacity of the quark-gluon gas. 
5.2 Heat Kernel on with Magnetic Field 

Let us now compute the heat kernel diagonal of the Laplacian with a constant magnetic field on R^. 
In the flat space limit the Laplacian on in polar coordinates with a constant magnetic field can be 
written as 

A = + , (5.20) 

where 

1 . „ r 
r 



Vi = cos ipdr sin ifd^- - sin LpH , 



V2 = sin ^ 5r "I — cos ipdui-\ — cos ipH . 
r 2 



(5.21) 
(5.22) 
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The heat kernel can be computed as follows. First, by observing that the covariant derivatives 
form the Heisenberg algebra 



[Vi,V2]=//, [Vi,//] - [V2,//] -0, (5.23) 
one can prove that the heat semigroup can be expressed in the form [2] 

exp(?A) - ^^-^ r dq exp(-i//cot(?//)M2 + (^,V>|, 



(5.24) 



where \q\^ = (q^)^ + (q^)^ and (q, V) = q^^i + q^^2- Further, it is not difficult to show that 

[exp {q, V> Six, x'yf"^ = 6{q) , (5.25) 
so that the integral over q becomes trivial and we immediately obtain the heat kernel diagonal [2] 

An sm(tH) 

Of course, since 'E? is non-compact, the spectrum of the Laplacian A is degenerate, that is, even 
if the eigenvalues are discrete their multiplicities are infinite. Therefore, the heat trace in the limit 
a — > 00 is infinite. However, the heat kernel diagonal is still well defined. Therefore, on the sphere S ^ 
in the limit as a — > 00 and fixed H = -X/a^ the heat kernel diagonal locally must have the following 
limit 

a->oo i 4ji sm{tH) 

Recall that // is a real anti-symmetric matrix with purely imaginary eigenvalues, so that this heat 
kernel diagonal is well defined. 

5.3 Heat Trace on 5 2 

5.3.1 Algebraic Approach 

Let {q^,q^,oS) be the canonical coordinates on the isometry group S0(3). Let C be the contour of 
integration in the complex plane of co defined by 

C = ^(C+ + C_) . (5.28) 

By using the isometry algebra and the representation of the Laplacian in terms of the twisted Lie 
derivatives one can show that the heat semigroup exp(fA) can be represented in form of an integral 
over the isometry group [7] 

C ^ ' R2 

x^^^^^exp[(^,X> + a;X3]. (5.29) 
0)12 
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where \q\^ = (q^)^ + (q^)^, {q, X> - q^-Ci + q^X.2- We will see later that one can take here either C+ or 
C_, which gives identical results. 

By using the isometry algebra one can show that [7] 



exp [{q, £) + ojL'i] 6{x, x') 



w/2 



^, sin [<x>/2] 

Substituting this into (5.29) we obtain the heat kernel diagonal in the form 



exp {ojG) 6(q) . 



U'^^^it) = exp 
47rt 



This can be written in the form 



ad(jj 



exp<- 



2 2 

a (jj 



ojll 



sin [a;/2] 



where with an integer n, is a function defined by 



f doj { 



6t» 



exp < 1- inci)!! 

' At 



w/2 



sin [<x>/2] ' 



Recall that G has imaginary half-integer eigenvalues. 

5.3.2 Properties of the Function ^(f; n) 
Let 

Jo* 



/do) { aP- , _1 
— : exp < 1- mo) 2 > 
V4^ 14? j 



w/2 



sin [(x>/2] 



so that 



By deforming the contour of integration one can show that 

n) = ^+{t, n) + i[R^{t; n) + «)] , 

where 



= (^)'''g(-l)^(i-)2;r)texp|-^} 
/?-(^;n) = (y) J](-lf>->2.fcexp 

A:=-l ^ 



(5.30) 



(5.31) 



(5.32) 



(5.33) 



(5.34) 



(5.35) 
(5.36) 

(5.37) 
(5.38) 
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Obviously, 
Therefore, 
and 



R^{t;n) = -R+{t;n). (5.39) 
^(f; n) = n) = n) , (5.40) 



^{t;n) = £^ exp{-^|-i^cos(«6>/2), (5.41) 

— DO 

where denotes the Cauchy principal value of the integral. This can also be written as 
nnn) = y (-1)«>-) f ^ e^pLUco^lnkA 



(a; + 27rtt)/2 

' cos(na;/2), (5.42) 
sin L<^/2J 



which is nothing but the sum over closed geodesies of 5^. 
We can also rewrite the function ^(f; n) in the form 



^ 2 

^(f; n) = -f exp ^^^^ , cos [nw ^^t/2\ , (5.43) 

J V4^ I 4 J sin[ojV?/2] ^ 

which allows the analytic continuation in the complex plane of t with a cut along the negative real 
axis (so that t = with \9\ < n), 

OO 

^(f; n) = exp l-^l "^^^^^ , cosh \nco V^/2l . (5.44) 

J V4^ I 4 J sinh[wV^/2] ^ 

Note also that the function *P(?; n) is an even function of n, that is, 

Y(f; «) - Y(f; -n) = ^(f; |?i|) , (5.45) 

so that it depends only on the absolute value \n\. In terms of the matrix G this means that *F(f; -2iG) 
depends only on the absolute value \G\ of the matrix G, which can be defined as the positive square 
root of the real symmetric matrix \G\ = V-G^. Recall that G has purely imaginary half-integer 
eigenvalues, therefore, \G\ has positive half-integer eigenvalues. 

We can also find the dual representation of the function ^'{t; n) as follows. As we have seen above 
one can use either contour, C+ or C for the calculation of the function ^{t; n). We use the series the 
series 

'^^'^ = -/oj V exp 1/ (z 4- ^ I ojI (5.46) 



sin[w/2] j-^ I \ 2 
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for Im a» > and 

coll 



sin [(x»/2] 

for Im a> < 0; they both converge absolutely. 
Now we notice that 

,2 



'^^_gexp|-/|/+^Ja;| (5.47) 



/d(jj ( o) ] toll 

—= exp <^ - — ^ . / sin[nw/2] = . (5.48) 
^|^t I 4fjsin[t^/2] 

C+~C- 

Indeed, the integrand is analytic on the real axis and, therefore, this integral vanishes. Therefore, we 
can add it to the original integral to get 

1 r du) { a>^l (jjjl 
'¥{t;n) = - —=tx^\-—\ . cos[?ifaj/2] 

c++c_ 

i C doj i (JJ^] ojjl . ^, , 
+ -^exp<^-— ^— — — sin[n|W2] 

c+-c_ 

I r doj ( ., , '^/2 

= - I — ^exp^j-— + iHcu/2 



V4^ I 4f J sin [aj/2] 

c+ 

I C dcj ( ., , <^/2 

— , exp< i\n\a>2> 

2j ^/4^t \ 4t /sin[o./2] 

C- 



Now, by using the expansion above we get 

,2 



Finally, by using the equation 

.2 



_ r ^/o. / oj\.( l + \n\\ \. 
——= exp < — — ±i\l + — - — \a)}iaj 

c± 

\2 



we obtain the spectral representation 



(5.49) 



(5.50) 



2.(/.l^ exp ^ (5.51) 



«F(, „) = . I 2 (/ . i^) exp |-. (/ . (5.52) 
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By using this equation we see that all functions ^{t, n) for different n can be reduced to just two 
functions, ^{t; 0) and 1), namely, for n > 



In) = »P(f; 0) - ? ^ (2/ + 1) exp |-f + ^) | ' 



^{t,2n + 1) = YCf; 1) -2?^ 



le' 



(5.53) 
(5.54) 



1=1 



Also, it is not difficult to see that the function *?(?; 0) can be expressed in terms of the function ^{t; 1) 
as follows 

^(f;0) = 24'|^;lj-^(f;l). (5.55) 

This property can be obtained directly from the integral representation (5.41). 

One can also show that the function ^'{t; n) has the following dual integral representation 



r / \n\\ \ I \n\\ 

'■¥{t;n) = t I dv tan(;ry)/ly + y lexp yl 



(5.56) 



where F is a contour that goes counterclockwise from (oo + /e) to (^ + is), then around the point ^ 
between and ^ on the real axis to the point (^ - is), and finally to (oo - ie). 

5.3.3 Heat Trace and the Spectrum 

By using the spectral representation of the function 1'(f; n) we can write the heat kernel diagonal of 
the Laplacian in the form 

^ oo 



+ |G| --+G' 



1=0 



where 



Ai{n) ^ l{l + n+l) + -, 
di{n) - 2l + n+ \ . 

Then the heat trace of the Laplace type operator L = -A + 2 is 



oo 

Tr exp(-fL) - tr ^ J, (2|G|) exp {--^ (2|G|) + o^q] 

1=0 ^ 



(5.57) 



(5.58) 
(5.59) 



(5.60) 



Recall that G is a generator of S 0(2) in some A'^-dimensional representation. We can diagonalize 
as follows 

^ /.'Ml .mN\ ..... 

G = diag (^jy,...,;— j , (5.61) 
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where mi, . . . ^m^ ai^e some integers. Since the matrix Q commutes with the matrix G, it can be 
diagonalized simultaneously with G, that is, 

Q ^ diag {Qu...,Qn) . (5.62) 
By using this decomposition the heat trace takes the form 

N CO 

Tr exp(-fL) ^ TjYj "^'^^""J^^ ^""P ("i^'-^) ' ^^"^^^ 

where 

- M\mj\) + a^Qj, j^l,...,N,. (5.64) 

Therefore, the eigenvalues of the operator L are -^^ij , j - I, ... ,N;l = 0,1,2, ... , with multiplicities 
di{\mj\). The smallest eigenvalue of the operator L is 

^min(Z.) - 4 ."li'l, l^oilmjl) + Qj] = mn + Qj\ ■ (5-65) 

a'^i<j<N^ ■'' \<]<N ya^ 2 J 

In terms of the function ^(f; n) the heat trace can be written as 

2 N 



Trexp(-?L)^-2]exp - -^-a^Qj 4'(-;my 



7=1 



(5.66) 



5.3.4 Asymptotics 

The asymptotics of the function 'F(f; n) as t ^ oo can be directly obtained from the eigenvalue repre- 
sentation; it is given by the lowest eigenvalue 

»P(f;n) ~(1 +|n|)fexp|- /^ "^J"'^ 1 + (3 + |»|) f exp |- /^ "^J"'^ } + ••• . (5.67) 



To compute the asymptotics as f — > we use the rescaled representation (5.43). This form is 
particularly useful to compute the short-time asymptotics as f ^ 0; we just expand the exponent in 
the powers of t and compute the Gaussian integrals over to. The asymptotics as f ^ are 

DO 

'F(f;?2) ~2Q(«)f\ (5.68) 

k=0 

where 

(5.69) 

Clin) ^ T^- '-r, (5.70) 



1, 






1 






12 ~ 


4' 




7 


«2 


«4 



Therefore, the heat kernel diagonal asymptotics as ? — > ai-e 

U'^'HO-i^nty'Y.^kQ , (5.73) 



k=0 

where 



1 1 1 2 

flo = 1, ai = -, ^2 ^ TT + tG . (5.74) 
i 15 D 

This coincides with the well-known general heat kernel coefficients (3.21). 

Moreover, this integral can be used to compute the asymptotics of the heat kernel diagonal as 

a — > oo. It amounts to just shifting the contour of integration in (5.31) by oj co + 2Gt/a^ to get 

oo 

,,diag,,x 1 ( t \ r doj ( aj^\ (jj^tl{2a) + Gtla^ 

U ^(t) = — exp — T- — — exp 

Ant \4a2 ; J V4^ \ 4 / sin [oj ^t/{2a) + Gt/a^] 

(5.75) 

Next, let G = r + a^H, where T = T(L) is a generator of the group S 0(2) in the representation T and 
H = -X(L)la^ is the generator of the group S 0(2) in the representation X taking values in the adjoint 
representation of the Lie algebra of the gauge group. Then 

DO 

,,diag,,x 1 I t \ r dw I o?\ tH + co^tl{2a) + tTla 
U^{t) - - — exp — -\ I — — exp -— ' 



4jTt V4a2 IJ^^ \ 4 / sin \tH + oj V?/(2a) + tTIa^] ' 

(5.76) 

Now, the asymptotics as a ^ oo can be computed by expanding the integrand in a Taylor series 
in inverse powers of a and computing the Gaussian integrals over u. This gives the correct leading 
asymptotics that coincides with the Euclidean space limit (with the magnetic field H) 

\imU^^Kt) = ^-^ (5.77) 
a^oo 4n sm {tH) 

Of course, such asymptotics describe the behavior of the heat kernel for large a and fixed t. The 
behavior of the heat kernel for large t is sensitive to the value of a, it is described by the lowest 
eigenvalue of the Laplacian, which is different on the sphere and on the Euclidean space R^. 

5.3.5 Spectral Approach 

As an alternative, let us compute the spectrum of the Laplacian directly. The Laplacian (4.44) is a 
differential operator acting on the space L?-{[0,an] x [0,2K],asm{r/a)drdip). We will find it use- 
ful to introduce a new variable x = cos(r/a). Then the Laplacian acts on the space L^{[-1, 1] x 
[0, 2n], a^dx dip) and takes the form 

A - ^ |(1 - x^)dl - 2xd, + [d^ + (1 - x)g]H . (5.78) 
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Let a,j6 be positive Hermitian operators defined by 



a = \d^\ , /3=\d^ + 2G\ , (5.79) 

and 

p{x) = (1 - xfil + xf . (5.80) 
We intertwine the Laplacian as follows 

A = p-'>^^p'l^ 

- -^(^a,/J + ^(a+yS+l)'-^+G'j , (5.81) 

where 

^Q-,/? = -(1 - -\fi-a-{a-r^ + 2)x] 5,, . (5.82) 

Then the operator A has the same spectrum and, therefore, the same heat trace as the original 
Laplacian. The operator Kafi is an elliptic self-adjoint operator on L?{[-1, 1] x [0,2k], p) with the 
weight function p{x). It is well-known [15] that the operator Kafi (with positive real a,yS > 0) has a 
discrete spectrum with simple eigenvalues 

AiiK^^p) = Kl + a + p + I) , (5.83) 

where / = 0, 1, . . . , with eigenfunctions proportional to the Jacobi polynomials p''"'^\x). 

This immediately enables us to compute the heat trace of a Laplace type operator L = -A -i- 2 as 
follows 

Tr exp(-?L) = fj^^' exp ||z + ^{a +/3 + 1)J +G^-^+ a^Q^j , 

(5.84) 

where the trace Tr in the last equation is the combined trace over the operator dip and the fiber trace 
over the matrix G. Since the eigenvalues of the operator are purely imaginary integers, im, with 
m € Z, this can finally be written as 

Tr exp(-fL) ^ tr exp \ - \ + a^Q 



OO CO 



/=0 m=-oa 



1 

/ -H - (|m| + \m - 2iG\ + 1) 



(5.85) 



where the trace tr is now over the fiber indices only. 

Now, one can show that this sum is equal to the previous result (5.57), (5.60), obtained by a 
completely different algebraic approach. 
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6 Effective Action 



6.1 Yang-Mills Effective Action 

The magnitude of the magnetic field in our scenario can be written as |;Fp - jK/a^, where K = 
' with p being the number of positive roots of the gauge algebra, and, therefore, the classical 
(Euclidean) action on the covariantly constant background S ' x S ' x 5 ^ is 



where cr = K/Se , 



16;rV— , 

xy 



a a 

y = 



(6.1) 



(6.2) 



ai aj 

and ai and a2 are the radii of the two circles. The total effective action is given by the sum of this 
classical part and the one-loop effective action (3.10). 

For the Yang-Mills ghost operator, the endomorphism Q and the generator G are 



tSghost ~ 



'-'ghost - ~^ ' 



(6.3) 



where X is the generator of the spin group Spin(2) taking values in the adjoint representation of the 
Lie group G. The volume of the manifold M = S^xS^xS^ is vol (M) = \67r'a^aYa2- By taking into 
account the contribution of the torus 5 ^ x S \ the heat trace of the Yang-Mills ghost operator is 



trAdCxp 



(6.4) 



For the Yang-Mills operator Lyect the endomorphism 2vect is given by 

a^Qyect = P + 2EX, (6.5) 
where P and £ are 4 x 4 real matrices of the form 

P"b = S"232b + S"^6'ib , E"h = S"2S3h - 5"i62b ■ (6.6) 
Notice that P is a 2-dimensional projection, that is, - P, and E-^ = -P. The generator Gyect is 

Gy,,t = E-X. (6.7) 

Therefore, 

^vect + ^^^Gvect = ■ (6.8) 

So, again, by adding the corresponding factors for the torus 5 ' x 5 ' we get the heat kernel diagonal 
of the Yang-Mills operator 

r , \ r , \ r , /I \i 

(6.9) 



0vect(O = 167r^fl:^fl:ifl;2n -T O -T trArfCxp 
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Thus, the total Yang-Mills heat trace is 



0m(O = l6K^a^aia2a{-^]n{-^]wYM(-^) , 



(6.10) 



where 



WYMit) 



r doj ( 



or] a>/2 



At I sin[tf;/2] 



xiVAd\ exp 



1 



t\--X^\-Xa) 



[tr Ti expC^o;) - 2] 



The trace over the vector representation can be easily computed. By using 

e.xp{EcjS) = 1- P + P cos oj + Esino), 

we get 



tr T, exp(£'(x>) = 2 + 2 cos oj . 



(6.11) 

(6.12) 
(6.13) 



Thus, the function WyAiit) is 



y/Ani 



expi--r- 



a>-] a>/2 



At I sin[tfj/2] 



2 cos ojtr Ad exp 



t\--X^\-Xu) 



Recall that the eigenvalues of the matrix X are imaginary half-integers 



Spec(X) = <|0,_^,/^,-/^,...,/^,-/^^,^ 



(6.14) 



(6.15) 



where p is the number of positive roots, r is the rank of the group, and ni, . . . ,np are some integers 
that we call monopole numbers. 

Therefore, the trace over the adjoint representation can be also computed to take more explicit 
form 



WyMit) = ^ exp 

7=1 



t(l+np 



r doj ( 

J vis '""l-- 

c 



uP' \ 0)12 



At I sin[w/2] 



4 cos ^njaj/2^ coso) 



-i-rexp - 



4/J ^/A^t 

c 



da> ( oj^] a>/2 

exp —. — ; — —2 cos a> , 



At I sin[oj/2] 



(6.16) 



which can be written in terms of the function *F(;; n) 



WYM(t) = 2 ^ exp 



til+np 



{Y(?; 2 + \nj\) + Y(f; 2 - \nj\)] + 2re"Mt; 2) . 



(6.17) 
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The ultraviolet (high-energy) properties of the effective action depend of the asymptotics of the 
function Wym as f ^ 0. By using the asymptotics of the function *P(f; n) we obtain 



WYM(t) ~ 2m + Cit + Cjt +■■■ , 



where 



C, 



7=1 



4 19 11 2 



(6.18) 

(6.19) 
(6.20) 



By using this result, we obtain easily 

T 1 / o A \ 

Pym = 167r — I C2 - za Ci + z a mj . 



xy 



(6.21) 



Now, by subtracting enough terms of the Taylor expansion at f = we obtain the renormalized 
heat trace 



4 

&^lt{a^t) = l6n'-[^l(xh)Q(y^t)WYM(t)-e-'"''RYMit;a^A)] , 



xy 



where 



(6.22) 
(6.23) 



RYM(t; a^^) = 2m + (Ci + Ima^Aj t + (C2 + Cio^A + ma^A^^ . 
Then the one-loop Yang-Mills elfective action is 

r(i)7M(-^'3'.j"'«^^z) = jl'-'^ ~ ^^^^1 -HzVm)log^ + (^YM[x,y;a^A,a^z)^ , (6.24) 



where 



dt 2 2 



(6.25) 



The convergence properties of this integral as f ^ 00 (infrared region, low energies) crucially 
depends on the asymptotics of the function WYM(t) as f ^ 00. Let 



Aoin) = 1 - 



1 _ M 
2 



( 2 for n = 0, 
- for = 1, 



-— for \n\>2. 
2 ' ' 



(6.26) 



Let 



Amin = min Ao{nj) 



(6.27) 
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and 

2 (l + |2 - \nj\ l) , (6.28) 



i<;<p 

''o(";'=''imn 

where the summation goes over all j such that Ao{nj) - /Imin- Then as f ^ oo 

WvMit) ~ 2?ce~'''™" + 6rfe"2' + • • • . (6.29) 

Thus we see that for \nj\ > 2 the lowest eigenvalue is negative, which leads to the exponential 
growth of the function WyMit) as f ^ oo and, as a result, to the infrared divergence of the integral for 
the zeta-function and the effective action and to the instability of the chiomomagnetic vacuum. The 
only stable configurations are those in which all monopole numbers are equal to or ± 1 . Recall that 
the flat space limit a ^ cxj is recovered slsX - -a^H — > oo, keeping H constant. That is, in the limit 
a ^ oo the matrix \X\ oo and, of course, the monopole numbers \nj\ oo (which can be interpreted 
physically as condensation of magnetic monopoles) leading to the instability. This is nothing but the 
instability of the Savvidy chromomagnetic vacuum in the flat Euclidean space. 

However, as we have shown above since the monopole numbers nj are not independent but are 
rather determined by the roots of the Lie algebra g it is impossible to have all monopole numbers 
bounded between - 1 and 1 . Therefore, the lowest eigenvalue is necessarily negative and is equal to 

'Imin = -2"max , (6.30) 

where ?imax - maxi<y<p \nj\ is the largest monopole number. Since n^^^ > 2, then the constant c is 
now 

C = (?lmax - 1)^, (6.31) 

where q is the multiplicity of the largest monopole number. 

Let us consider the group SU{N) for simplicity. Then the monopole numbers are riij = ki - kj, 
where ki are some integers whose sum is equal to zero. For N = 2we have just two non-zero integers, 
ki and k2, whose sum is equal to zero. Then there is one non-zero monopole number 

ni = 2ki . (6.32) 

Therefore, 

"max ^ 2\ki \ > 2 and ^min = -1^1 1 < -1- (6.33) 

For N = 3 there ai^e thi^ee non-zero integers ki,k2, k^ whose sum is equal to zero. Then there ai^e 
three non-zero monopole numbers 

ni = ki - k2 , «2 - 2^1 + k2 , nj, = ki + 2^2 ■ (6.34) 

Then, if ki = k2, 

«max-3|^l|>3. (6.35) 

If ki > k2 > 0, then 

"max = 2^1 -I- ^2 > 3 . (6.36) 
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Ifk2 <0<2|yt2l <ki,±en 

"max = 2^1 - 1^2! > 3 . (6.37) 
lfk2<0 <ki and \k2\ < 2ki < 4\k2\, then 

"max - ^1 + \k2\ > 3 . (6.38) 

Finally, iflfk2<0< ki and 2ki < \k2\, then 

«max-2fe|-^i >3. (6.39) 
In this case ^ 

^min<-2- (6-40) 

More generally, one can show that for the group S U (N) we have n^nax ^ 2 for even A'^ and 
"max ^ 3 for odd A'^. The equality is achieved by choosing the non-zero integers ki,. . . ,kj^, (whose 
sum is equal to zero) such that there are as many as possible (monopole-antimonopole) pairs (+1, -1). 
Then for even N, the largest difference \ki - kj\ is equal to 2. For odd A'^, we will have as many pairs 
(+1, -1) as possible and a triple (2, -1, -1). Then the largest difference \ki - kj\ is equal to 3. Thus, 
the bottom eigenvalue for the group SU{N) is less or equal to (-1) for even N and less or equal to 
(-3/2) for odd A^. 

Let us compute this function for the group SU{2N) so that ni - 4N^ - I, r - 2N - 1 and 
p = N{2N - 1), in the simple case when the numbers kj are chosen as A^ pairs (+1,-1), that is, 
ki = ■ ■ ■ = k^ = I and ki\f+i - • • • - k2N = -1- Then there are only A'^^ nonzero monopole numbers 
which are equal to 2 and the rest are equal to zero. Then 

^FM ^^^(^) = ^N^e^"^ {^fCf; 4) + ^(f; 0)1 + 2i2N^ - l)e"'^Wit; 2) (6.41) 
= t [2N^e' + 6N^e~' + 6{2N^ - l)e-^'] 
+ [4N^e^"'^ + 2{2N'^ - l)e'l'^] »P(?;4) . 

The asymptotics of this function as f — > cxd are 

Wl^^^^\t) = 2N^te' + 0{te-') . (6.42) 

A few remarks are in order. First of all, the effective action does not depend on the parameter A; 
it can be chosen arbitrarily, for example, A = l/a^. Second, for a sufficiently large Rez the function 
^YM{x,y;a^A,a^z), as a function of z, is analytic. Since the function WyMiO grows exponentially as 
f ^ 00, the function <^ym has a branching singularity at 

ZO = 5- = T-Tr"max ■ (6.43) 

2a^ 
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6.2 Effective Action of Matter Fields 

For the scalar operator the potential term is Qq, which is determined by the second derivative of the 
scalar potential V{(f>), and the generator Go is 



(6.44) 



where Xq is the generator of the spin group Spin(2) taking values in the representation Wq of the gauge 
group. Therefore, the heat trace of the operator Lq is 



0o(O = l6Tr'a^a\a20. 1 I O I I tr exp 



--{Xl + a'Q,-]^ 



^(^;-2/Xo) 



For the spinor operator the generator Gspin is now 

1 



G 



spin ~ 2 ^ -^spin i 



(6.45) 



(6.46) 



where y - y^.j/^, with ja being the Dirac matrices, so that is the generator of the spinor represen- 
tation; of course, - -I, and Xgpin = X^^\a(^M) is the generator of the spin group Spin(2) taking 
values in the representation Wspin of the gauge group. The endomorphism 2spin is 



1 



1 



Gspin = + ^T^spin + M 



Then 



G^pi„ -Ha^ Gspin 



7l + 4in + «'^'- 



Therefore, the heat trace for the operator Lgpin is 



(6.47) 
(6.48) 



0spin(O = l6Tr'cP'a\a20. 



t 

(a J 



— Ifl 



^^rw.p.„trr,p,„exp 



x'F(^;/(2Xspin-r)) ■ 



Thus the total heat trace for the matter fields reads 



0mat(O = 167r3a2^ia2n[4]n[4]w^mat(4) ' 



(6.49) 



(6.50) 



where 



r dw ( 



CO 



oj/2 



At I sin[if;/2] 



xariv„exp 



-t\Xl + a^Qo--\-Xoaj 



^Kin + «'^') + (^r-^spin)^ }. (6.51) 
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To proceed further, we assume that the matrix M does not transform with respect to the spinor 
representation (it does not have spinor indices, but only the group indices). Then by using 



exp 



= cos(tL>/2)I + y sin(6L)/2) , 



the trace over spinor representation is easily computed 



Thus, we get 



Wmat(0 



trj- exp 

^ spin ^ 



/do) { 



1 



CO 



4cos(aj/2) . 



cj/2 



(6.52) 



(6.53) 



4t I sin[a>/2] 



-t\Xl+a^Qo--\-Xoaj 



x|triv„exp 

-4tr w,p.„ exp [-t (x^pj^ + a^M^) - X^pinO;] cos(a;/2) 
The eigenvalues of the matrices Xq and Xgpin must be imaginary half-integers 

Spec(Xo) = {/^,...,/^,}, 



Spec(Xspin) = 



(6.54) 



(6.55) 



(6.56) 



2 2 

where m, and kj are some integers determined by the weights of the representations reaUzed by spinors 
and scalars. 

For simplicity, we assume that the symmetry is not broken so that all scalar fields have the same 
mass niQ and all spinor fields have mass m^pin- Then the potential terms Qq and M are proportional to 
the identity, that is. 



Qo = /MqI , M = m 



spin J" 



(6.57) 



In this case 



WmAt) = exp 



2 7 

— - a I ? 



No 



Y^e""V^^>(t;mj) 



-2e-'"'<'^ £ e"^'^ [^(t; 1 - kj) + T(?; 1 + kj)] . 



(6.58) 



The leading asymptotics of this function as ? ^ oo is 



Wr^M ~ rote-"'''^ - 2 J] \kj\ te~'"''<'^^ , 



(6.59) 
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where rg is the number of zeros among the scalar weights mj. 

This indicates infrared instability for massless fields, when mo - 'Mspin = 0. It is interesting to 
note that the infrared instability of massless spinor fields is not caused by the zero numbers kf, in fact, 
that contribution is proportional to ^{t; 1) and has a nice exponentially decreasing behavior at ? — > oo. 
This instability is intrinsic and cannot be avoided if there are non-zero numbers kj. This simply means 
that massless Dirac operator with a constant non-zero magnetic field always has zero modes onS^. 
We can even compute the multiplicity of the zero eigenvalue, that is, the dimension of the kernel of 
massless Dirac operator. 



dim Ker L, 



spin 



-2£|^;| = 4tr 



Wspinl^spinl 



(6.60) 



Let us compute the contribution of spinors in the fundamental representation of the group S U (2N) 
when the numbers ki, . . ., k2N ai"e chosen as ki = ■ ■ ■ = kN = ^ and k^+i - • • • = k2N = -1, so that 



A^,pi„ ^ 2N. Then 



spin ^ ' ^ 



-ANte 



4 - «^'"spin 1 1 



- SA^exp 



m?;2) + 'F(f;0)] 



4 - ^^'"spin I f 



*P(?;2). 



(6.61) 



Now, the heat trace for matter fields can be renormalized in the same fashion as we have done for 
Yang-Mills fields. Then we get the effective action for matter fields exactly in the same form as for 
the Yang-Mills fields with the function <^max{x,y,a^A,a^z) that is expressed in terms of the function 
W^mat(0 exactly in the same way as for the Yang-Mills fields. 



7 Thermodynamics of Yang-Mills Theory 

In this section we investigate the entropy and the heat capacity of the gluon gas. The temperature is 
related to the radius of the first circle by T = ^^^^ ■ The volume of the space is expressed in terms of 
the radius of the second circle and the radius of the sphere V - SK'^a2a^ ■ 

For a canonical statistical ensemble with fixed T and V the free energy F = E - TS is a. function 
of r, V defined by 



Then the enti'opy is defined by 



1 r 

F - rr = 

In ai 



d . d T 
S = -—F = a^- 



dT 



and the heat capacity at constant volume is 



(7.1) 



(7.2) 



^ T—S = -a] 
dT ' 



52 F ^ 5 F^l 

fll T h 2 



(7.3) 
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By using our results for the effective action we obtain 

ai 2a^ [ A \ai 02 J) 

where the function O is given by the sum of the corresponding functions for the Yang-Mills fields and 
the matter fields, O = (^ym + ^mat, and 



1 

167r3 aia2 



(J3ym + J3mat) . (7.5) 



It is easy to see that neither the classical part nor the coefficient b depend on the temperature. 
Therefore, the entropy and the heat capacity do not depend neither on the classical term nor on the 
renormalization parameter /u. Recall also that the parameter A is completely arbitrary; so we can set 
it equal to A = l/a^. Thus the entropy and the heat capacity of the quark-gluon gas are given by the 
derivatives of the functions O 

S ^ ^^Ax,y,l,ah), (7.6) 
2y 

Cv - ~0,Ax,y, ha^z) - ~<^xAx,y, hah). (7.7) 
2y 2a\ 

The entropy and the heat capacity per unit volume are 

^ = -^^Ax,y,hah), (7.8) 

V = T7-^^--(^, y, 1, ah) . (7.9) 
Now, by differentiating the function O with respect to x we get 

—e~'"''a' [x^t) a (y^i) W{t) , (7.10) 
je~'"'' [n' [x^t) + 2}? to." (x^t)] Q. (y'^t) W{t) 

= 2^f f.-"'--|a'<o.2,a"<o|a(5,),v(l). 



where 



(7.11) 



W{t)^WYM{t) + W^,,{t). (7.12) 



Now by using the asymptotics of the function Q. we easily see that the integrals (7.10) and (7.1 1) 
converge at f ^ 0. Moreover, for large z these integrals also converge as f ^ 00. Since the function 
Wym increases exponentially at infinity, the function <1) has a singularity at a finite positive value of z- 
However, because of the asymptotics of the function Q.'{t)+2t£l"{t) as f ^ 00, (5.19), we immediately 
see that the integral (7.11) for the heat capacity may converge even in the infrared limit z ^ 0, due to 
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the presence of an extra exponential factor e~' . Namely, for the groups SU{2N) there are monopole 
configurations such that the function WyMit/x^) increases as e'^^^ , and therefore, the integral would 
converge for sufficiently large x. 

That is why we investigate this case in more detail. We set z = (and y = 0, for simplicity) to get 

-[a'{t) + ita''{t)]wy-Y (7.13) 

We consider the group SU(2N) when the function W is given by (6.42) and (6.61) for the spinors. 
We decompose it according to 

where = for massive spinors and * = 1 for massless spinors, when mspin = 0, and the function V 
is exponentially decreasing as f ^ oo. 

We now split the integral (7.13) into three parts accordingly 

0.„(^, 0; 1,0) = /i + /2 - SNsvi , (7.15) 

where 

r°° dt 

vi = — {a' (t) + 2ta" (t)} , (7.16) 

Jo ? 

and 

/i = 4^2 / _^t/x' j^/ (-J) ^ 2ta" it)] , (7.17) 
Jo ^ 

-|Q'(0 + 2fn"(0|v(— j . (7.18) 

Recall that x = a/ai = 2naT; so ;c ^ oo is the high-temperature limit and x — » is the limit of zero 
temperature. 

We consider the high-temperature Umit first. As ;c ^ cx) the second integral is 

I2 ~ 4(4N^ - 1)V2X^ , (7.19) 

where 

V2= \ -T {n' (0 + 2ta" (01 . (7.20) 

Jo f 

For the first integral in the limit x ^ 00 we get 



Thus, as r ^ 00 we obtain 



/i ~ 4n\i . (7.21) 



r 

y ~ 27r2(4Af2 - 1)V2T^ + ■■■ , (7.22) 



which is similar- to the black body photon gas. 
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The first integral exliibits mucli more interesting behavior due to the presence of the growing 
exponential factor e'l'^^ . It converges at f ^ for any x. However, at ? ^ oo its convergence depends 
on the critical temperature. Recall that the other part of the integrand has an exponential factor e~' . 
Therefore, the integral converges in the infrared domain ? — > oo if ;c > 1 and diverges otherwise. This 
defines the critical temperature 

Tc = :^. (7.23) 
2na 

At the critical temperature T k the second integral is simply 

/2 = 2 -j[n'{t) + 2tn" (?)} V (0 . (7.24) 
Jo f 

However, the first integral is nonanalytic neai^ T^, namely, by using the asymptotics (5.19), we 
get as X ^ 1^ 

/i~8A^^fl-ir'~^(x-l)-^/^ (7.25) 



x^J V2 

Therefore, the heat capacity near the critical temperature, T ^ , is 

Cv _jf_lTj-TX^''^ 
'^~47rV2a3\ T, ) 

This indicates the second-order phase transition with the critical exponent a - 3/2. 



(7.26) 



8 Conclusion 

The primary goal of this paper was to study of the Yang-Mills vacuum in the low-energy (long- 
distance) limit. The Savvidy model of such a vacuum with constant chromomagnetic field in Minkowski 
spacetime suffers from a well-known instability, which exhibits itself in negative eigenvalues of the 
gluon operator. We noticed that a positive spatial curvature of the spacetime manifold acts as an ef- 
fective mass term and, therefore, can stabilize the Savvidy vacuum. That is why we considered the 
case of a compact spacetime manifold of the form S ' x 5 ^ x S ^ with a covariantly constant chro- 
momagnetic Yang-Mills field on the sphere S^. On the sphere, such a configuration is of monopole 
type parametrized by a collection of half-integers (monopole numbers). Such a configuration has a 
well-defined (Savvidy type) flat space limit with constant chromomagnetic field when the radius a of 
the 2-sphere S^, as well as the monopole numbers, nj, go to infinity, nj,a — > oo, such that the ratio 
rijla^ (which defines the magnetic field) remains constant. This limit can be interpreted physically as 
the condensation of monopoles. 

We computed exactly the spectra and the trace of the heat kernels of all relevant operators, which 
enabled us to compute exactly the one-loop effective action. We have found that the gluon operator 
does not have negative eigenvalues only when the monopole numbers are between -1 and 1. That is, 
any monopole number Uj, with \nj\ > 2, leads to an instability of the chromomagnetic vacuum. This 
confirms once again that the flat space limit with constant chromomagnetic field is unstable since it is 
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created by infinitely large monopole numbers, formally \nj\ oo. We showed that for any compact 
simple gauge group there are always monopole numbers with absolute value greater or equal to 2, 
which means that there is no stable constant chromomagnetic configuration also in curved space (at 
least in the model 5 ^ x 5 ' x S^). 

We also studied the thermal properties of Yang-Mills theory, in particular, we computed the 
entropy and the heat capacity of the quark-gluon gas. We have found that the heat capacity is well 
defined even in the infrared limit and computed the high-temperature asymptotics of the heat capacity. 
Moreover, in a particular model S U (IN) we found that the heat capacity has a typical branching 
singularity ~ (T - TcY^^^ at a finite critical temperature Tc = ll{2na) indicating the second-order 
phase transition. 

We conclude that to stabilize the chromomagnetic vacuum at lower energies one should consider 
non-constant magnetic fields on non-compact spaces. Constant magnetic fields on compact symmetric 
spaces are too rigid, they are completely determined by the spin connection and are of the same order 
as the space curvature. This makes it impossible that the gluon operator with the potential term 
R"b - 1T"b is strictly positive. What one needs is a large Ricci tensor and a small independent 
magnetic field to make this work. 

It is also interesting to study the Yang-Mills vacuum on the Einstein model 5^x5^. We intend to 
caiTy out such a study in future work. Another interesting question to pose is whether the Yang-Mills 
vacuum is stable when gravity is treated as a dynamical field. The technical calculations are not that 
difficult, but then, of course, we immediately face the non-renormalizability of general relativity, so 
we have to consider higher-derivative gravity instead. 

It is hard to imagine that this model can be directly relevant in hadron physics in the study of the 
confinement because of the completely different energy scales dictated by the gravitational constant 
and the cosmological constant. However, it can be relevant in the study of the structure of the quark- 
gluon plasma in the early Universe. 
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